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Abstract 

The first observation of the B® — > rfrf decay is reported. The study is based 
on a sample of proton-proton collisions corresponding to 3.0 fb -1 of integrated 
luminosity collected with the LHCb detector. The significance of the signal is 6.4 
standard deviations. The branching fraction is measured to be [3.31 ± 0.64 (stat) ± 
0.28 (syst) ± 0.12 (norm)] x 10~ 5 , where the third uncertainty comes from the B± —> 
rjK ^ branching fraction that is used as a normalisation. In addition, the charge 
asymmetries of B± —> rj'K± and B± —> which are control channels, are 

measured to be (—0.2 ± 1.3)% and (+1.7 ± 1.3)%, respectively. All results are 
consistent with theoretical expectations. 
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Hadronic decays of beauty hadrons into final states without charm quarks (charmless 
decays) are suppressed in the Standard Model of elementary particles. They proceed 
predominantly through b —>■ u transitions, mediated by the emission of a virtual W ± boson, 
and b —>• s transitions, mediated by the exchange of a virtual W ± boson and a virtual 
quark. The respective “tree” and “penguin” amplitudes are of similar size, allowing for 
possible large quantum interference effects measurable as charge-parity (CP) violating 
asymmetries. New particles not described in the Standard Model may contribute with 
additional amplitudes, and therefore affect both the decay rates and the CP asymmetries IT). 
The B ± — y rj'K ± and B° —>■ rj'K 0 decays, 0 first observed by the CLEO collaboration |2|, 
have some of the largest branching fractions among all charmless hadronic H-meson 
decays |3|. Studies of such decays, conducted so far mostly at e + e~ colliders operating 
at the T(45) resonance, provide accurate measurements of integrated (4j 51 and time- 
dependent j6j 7| CP- violating asymmetries in charmless hadronic B ± and B° meson decays, 
and are useful to look for deviations from Standard Model predictions. 

Charmless hadronic B® decays are poorly known, in particular decays to a pair of 
unflavoured neutral mesons |8j|9), but have been extensively studied in the framework of 

perturbative QCD 


QCD factorisation 10-12 


13 , soft collinear effective theory 14 , and 


flavour SU(3) symmetry 15 . The decay B s °4 rfrf is expected to have a relatively large 
branching fraction, similar to that of its SU(3) counterpart B —y rfK ; predictions range 
between 14 x 10~ 6 and 50 x 10~ 6 , and have large uncertainties 10 15 . The rfrf final state 


is a pure CP eigenstate. Decays to this final state of B° s and B° s mesons flavour tagged 
at production may therefore be used to investigate time-dependent CP asymmetries in a 
complementary way to the measurements in B® —> (jxj) 16 , but without the need for an 
angular analysis. 

In this Letter we present the first observation of the B° s —» rfrf decay. Its branching 
fraction is measured using the known B ± —y r)'K ± and B ± —> tf>K ± decays as calibration 
channels. The CP asymmetries of the calibration channels are also measured, relatively 
to the B ± —y J/xf> K ± channel. All these measurements use proton-proton (pp) collisions 


corresponding to 3.0 fb -1 of integrated luminosity, of which 1.0 (2.0) ft> was collected in 
2011 (2012) at a centre-of-mass energy of 7 (8) TeV with the LHCb detector. 

The LHCb detector 17 is a single-arm forward spectrometer at the LHC covering the 


pseudorapidity range 2 < 77 < 5, designed for the study of particles containing b or c quarks. 
The detector includes a high-precision tracking system, two ring-imaging Cherenkov 
detectors used to distinguish different types of charged hadrons, a calorimeter system 
consisting of scintillating-pad and preshower detectors, an electromagnetic calorimeter and 
a hadronic calorimeter, and a muon system. The trigger (18] consists of a hardware stage, 
based on information from the calorimeter and muon systems, followed by a software stage, 
which applies event reconstruction using information from all the detector sub-systems. 

Signal -Bj?—)• 7 / 77 ', B ± — > rfK , and B ± —> (f)K ± candidates are reconstructed through 
the decays rf — > 7r + 7r _ 7 and (jr —> K + K~. Selection requirements are chosen to be as 


1 Charge conjugation of neutral Z!jL mesons is implied throughout this Letter. The notations rf and </> 
refer to the if (958) and 0(1020) mesons, respectively. 
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similar as possible for the three channels and are optimised for B®—} rfrf , maximising the 
figure of merit e/(a/ 2 + '/B) [19], where e is the efficiency for selecting simulated signal 
events, B is the number of background events in the signal region estimated from the mass 
sidebands, and a = 5 is the target significance of the possible signal. The requirements 
on the 0 meson and on the charged kaon associated with the rf or 0 resonance in the 
candidate B ± decays, which is referred to as the bachelor kaon, are chosen to minimise 
the relative statistical uncertainty on the B ± —> 0/1 : ± signal yield. 

Charged particles are required to be inconsistent with originating from a primary 
pp interaction vertex (PV) and to have a transverse momentum (pp) with respect to 
the beam line in excess of 0.25GeV/c, while bachelor kaons must have p r > 1.2GeV/c. 
Particle identification algorithms are applied to distinguish kaons from pions, and photons 
from electrons 20 . Photons are required to have p p > 0.5 GeV/c. The intermediate rf 
( 0 ) resonances must have pp > 1.5(0.5) GeV/c and momentum p > 4 GeV/c. The 7r + 7r~ 
invariant mass in candidate rf decays must exceed 0.56 GeV/c 2 . 5-meson candidates 
must have pr > 4 GeV/c. Topological variables are used to isolate the signal, such as the 
angle between the reconstructed B momentum and the vector pointing from the PV to 
the B decay vertex (required to be smaller than 10 mrad), and the distance of closest 
approach to the PV of the B trajectory (required to be less than 0.04 mm). Reconstructed 
invariant masses of the B®, B ± , rf , and 0 candidates are required to be in the ranges 
5000 < m B o < 5600 MeV/c 2 , 5000 < m B ± < 5500MeV/c 2 , 880 < rrty < 1040MeV/c 2 and 
1000 < r?7.0 < 1050 MeV/c 2 , respectively. Only candidates with a well reconstructed B 
decay vertex are retained; in the events with multiple candidates (< 5%), the candidate 
with the smallest vertex y 2 is kept. 

The Bf —> rfrf signal yield is determined from a multidimensional unbinned ex¬ 


tended maximum likelihood fit to the sample of Bf —> 
dates, using the combined y/s = 7 and 8TeV data sets. 


and B ± —» rfK ± candi- 
The likelihood is written as 


rfrf 


C = exp(— V . Nj) 


n 


N 


V N i p, i 


where Nj is the yield of fit component j (signal or 


backgrounds), Pf is the probability of event i for component j, and N is the total number 
of events. The probabilities Pfj are expressed as products of probability density functions 
(PDF) for the invariant masses used as observables in the fit: the rfrf invariant mass 
( m r?v), ^ ie fw° randomly ordered 7r + 7r _ 7 invariant masses (m( 7r ^- 7 ) 1 and m^ ny ) 2 ) of the 
Bg ~* r f r f candidates, and the r) , K ± and 7r + 7r _ 7 invariant masses (rn ri >k and m W7r7 ) of the 
B ± —> tfK ± candidates. In the reconstruction of the 7r + 7r _ 7 candidates, the known rf 
mass [3] is applied as a constraint to calculate the and m^K variables. 

The B ± —>■ rfK ± sample is described with three components: the signal, and two 
combinatorial background components with and without an rf resonance in the decay 
chain. The Bf —> rfrf sample is modelled with seven components, three of which are 
significant: the signal, the combinatorial background and partially reconstructed 6-hadron 
decays without rf resonances in the final state. The remaining backgrounds, for which the 
event yields are found to be consistent with zero, consist of two combinatorial and two 
partially reconstructed components, each involving only one resonant 7r + 7r“7 candidate. 

All the PDFs that peak at the B or rf mass are modelled by Crystal Ball (CB) 
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functions 21 modified such that both the high- and low-mass tails follow power laws and 


account for non Gaussian reconstruction effects. The parameters used for the description 
of m v 'K for the B ± signal are free in the fit, while all the parameters of for the 
B° s signal PDF are determined from simulation, except the CB width. The ratio of the 
CB widths in m v ’K and m,y v > is fixed to that measured in simulation. The partially 
reconstructed background is described with an ARGUS function 1221 convolved with a 
Gaussian resolution function of the same width as the corresponding signal PDF, while 
the combinatorial background is modelled with a linear function. A common CB function 
is used for modelling all the rf resonances, with mean and width free in the fit, while tail 
parameters are determined from simulation. The mass distribution of r/' candidates from 
random combinations is modelled with an empirical quadratic function. 

We observe 36.4±7.8 (stat) ±1.6 (syst) B ° —>■ rfrf decays corresponding to a significance 
of 6.4 standard deviations, including both statistical and systematic uncertainties, discussed 


later. The significance is computed using Wilks’ theorem 23 , and is scaled by the ratio of 


the statistical over the total uncertainties. The measured B ± —> r]'K ± yield is 8672 ± 114, 
where the uncertainty is statistical only. Mass distributions are shown in Figs, [l] and [2j 
with the fit results overlaid. 

To measure the ratio of the branching fractions *6(5°—>• rf rf) / B(B ± —> rfK ± ), the fit 
is repeated taking into account the different reconstruction efficiencies in the 7 and 8 TeV 
data sets. The four background components with yields consistent with zero are neglected 
in this case. The common parameters between the 7 and 8 TeV data sets are the shape 
parameters and the ratio of branching fractions. The ratio of branching fractions is related 
to the ratio of yields according to 


B(B° a ^rfrf) f d 1 

B(B ± ^ifK ± ) f s B(r)' —>• 7 r+ 7 r _ 7 ) 
N t (B° s -> rfrf) w £ Z (R ± -> rfK±) 

X rfK±) X e^B^rfrf) ’ 


( 1 ) 


where the subscript l indicates the 7 TeV or 8 TeV data set, the ratio of probabilities for a 
b quark to produce a Bf or B ± meson is f s /fd = 0.259 ± 0.015 [24], and the branching 
fraction of the rf decay is B{rf —>■ 7 r + 7 r _ 7 ) = 0.291 ± 0.005 |3|. The ratio of efficiencies for 
reconstructing the normalisation and signal decay channels £i(B ± —, y rfK ± )/ei(B < f —)■ r]'r]') 
is determined from control samples (particle identification, photon reconstruction and 
hardware trigger on the signal) and simulation to be 8.46 ± 0.35 for the 7 TeV and 
7.85 ± 0.26 for the 8 TeV data sets, including all experimental systematic uncertainties. 
The largest uncertainty in the determination of the efficiency ratio comes from the photon 
reconstruction efficiency. This efficiency is measured using B ± —y J/if K*± decays followed 
by J/if —y /i + /± and K M± —y K ± tc° K ± 77 1 25 1 , and a crosscheck is provided by the 


measurement of the ratio of branching fractions of the B ± —» rf K ± and B ± —y cf)K ± control 
channels. The result of Eq. Q is then 


B(B° S ^ rfrf) 
B(B±^rfK±) 


= 0.47 ± 0.09 (stat) ± 0.04 (syst). 
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Figure 1: Distributions of the (top) rjrj, (bottom) (77777)1 and (77777)2 invariant masses for the 
Bg —> rjrj candidates with fit results overlaid. The (77777)1 and (77777)2 mass distributions are 
shown for the candidates with an 7/7/ invariant mass within three standard deviations of the B® 
mass. The components are the following: (dashed red curves) £7 —> rjrj signal, (long-dashed blue 
curves) combinatorial background without an 7/ resonance in the final state, (dot-long-dashed 
black curves) partially reconstructed background without an rj resonance, (short-dashed red, 
short-dashed blue curves) combinatorial background with one rj resonance, and (dot-dashed 
green, dot-dashed black curves) partially reconstructed background with one rj resonance. The 
total fit function is shown as the solid blue curves. 


Contributions to the systematic uncertainties are summed in quadrature leading to a 
total systematic uncertainty on the B ° —> ?/?/ signal yield (ratio of branching fractions) of 
1.6 (0.041). The uncertainties on f s /fd, B (7/ — >• 7 r + 7 r _ 7 ), ei(B ± — > 7/ K^) / £i(B° s — >■ 7/7/), 
and on the values of fit model parameters fixed from simulation, lead to a systematic 
uncertainty of 0.7 (0.038), while a variation of the PDF models leads to an uncertainty 
of 1.4 (0.007). The fit bias, evaluated in simulation, is consistent with zero, and its 
statistical uncertainty of 0.4 (0.005) is applied as a systematic uncertainty. Finally, an 
uncertainty of 0.014 is assigned to account for the neglected background components in 
the branching-fraction fit. 

Using the known value B{B ± —>■ rj'K ± ) = (7.06±0.25) x 10~ 5 |3|, the branching fraction 
is measured to be B(B® —* rj 7/) = [3.31 ± 0.64 (stat) ± 0.28 (syst) ± 0.12 (norm)] x 10 -5 , 
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Figure 2: Distributions of the (top) rfK ^ and (bottom) 7r + 7r - 7 invariant masses for the B± —> 
rfK± candidates with fit results overlaid. The components are the following: (dashed red curves) 
B± —> rf signal, (long-dashed blue curves) combinatorial background without an rf resonance 

in the final state, and (dotted blue curves) combinatorial background with an rf resonance. The 
total fit function is shown as the solid blue curves. 


where the third uncertainty comes from the B ± —> rf K ± branching fraction. 

The B ± —» rfK ± and B ± —> <pK ± charge asymmetries, A CP = (T~ — T + )/(r _ + T + ), 
where is the partial width of the B ± meson, are determined using the strategy adopted 
in Ref. |26|. For these measurements, we consider either events triggered on signal 
candidates (TOS events) or events triggered at the hardware stage independently of the 
signal candidate (non-TOS events). The raw asymmetries are obtained from a fit to 
the positively and negatively charged candidates, and each subsample is further split into 
TOS and non-TOS events, to account for trigger-dependent detection asymmetries 26 


For each channel a two-dimensional fit of the mass distributions of the B candidate and 
its neutral daughter is performed, and the four samples are fitted simultaneously. 

The model in the B ± —>■ r) , K ± fit is the same as that used in the simultaneous B° s —> rfrf 
and B ± —>■ rfK ± fit; the B± —> (j)K ± model mirrors that of Ref. 26 , except for the ul^k 
signal PDF, which is described by the sum of a CB function and a Gaussian function. The 
shape parameters used to describe the peaking and partially reconstructed component 
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Figure 3: Distributions of the (top) (f>K± and (bottom) K + K~ invariant masses for the B± —> 
4>K± candidates with fit results overlaid The components are the following: (dashed red curves) 
B± —> (f>K± signal, (long-dashed red curves) non-resonant B± —> I\ + K~ background, partially 

reconstructed 6-hadron background (long-dashed blue curves) with or (dot-dashed blue curves) 
without a 4> resonance in the decay chain. The combinatorial background with or without a cj) 
resonance is too small to be visible. The total fit function is shown as the solid blue curves. 


PDFs are shared, while independent parameters are used for the combinatorial component 
PDFs in the trigger sub-samples. In both fits, each component has four yields parametrised 
as N± = N k ( 1 =F where N k and A? p k are the yield and the raw asymmetry in 

each trigger category k, respectively. 

The observed B ± —» 7]'K ± [B ± » 4>K ± ) raw asymmetries and their statistical uncer¬ 
tainties are -0.019 ± 0.014 (+0.003 ± 0.014) and -0.027 ± 0.020 (-0.011 ± 0.018) for 
the TOS and non-TOS categories, respectively, the fraction of event in the TOS category 
being 63.8% (60.1%). The fitted mass spectra for the B ± —> <f)K ± candidates are shown in 

Fig- H} 

In order to determine the CP asymmetry, the raw asymmetry is corrected for the 
production asymmetry in pp collisions, ^4p, and for the bachelor K ± detection asymmetry 
due to interactions with the detector matter, AB,k- Under the assumption that these 
asymmetries are small, A ^ k is related to the CP asymmetry A CP as A^ p w k = A CP + 
^4d,A: + Ap. Because the CP -violating asymmetry in is known precisely |3|, 
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the raw asymmetry in this decay is used to determine the sum of the detection and 
production asymmetries. The raw asymmetry of the B ± —y J/iliK ± decay is measured in 
a simultaneous fit- to the mj/^K ± distributions of the positively and negatively charged 
candidates selected with similar criteria as for the signal modes. Independent fits are 
performed for events belonging to each trigger category. The fit model consists of a peaking 
signal component, described with the sum of a CB function and a Gaussian function, and 
a linear combinatorial background component. 

The raw B ± —> J/ipK ± asymmetries, —0.020 ± 0.004 (stat) and —0.011 ± 0.003 (stat) 
for the TOS and non-TOS categories, respectively, are subtracted from the B ± —)■ rj'K ± 
or B ± —» (f)K ± raw asymmetries for each trigger category. The weighted average of 
these asymmetry differences, AA CP = A CP — A CP (B ± —> J/ijjK ± ), is computed with 
weights given by the fractions of signal events in each of the two categories. The resulting 
asymmetry differences are AA CP (B ± —> rj'K ± ) = —0.005 ± 0.012 (stat) and A A CP (B ± —> 
0 K±) = +0.014 ± 0.011 (stat). 


Three significant sources of systematic uncertainties are identified. The first accounts for 
the effect of the mass shape modelling, leading to an uncertainty on AA C ' P of 0.021 x 10~ 2 
(0.20 x 10“ 2 ) for the B ± —> \rj’K± (B ± —> 0A' P ) channel. To account for the different- 
kinematic properties between the signal channels and the B ± —s- J/ip it ± channel, AA CP is 
measured in three independent subsamples selected according to the transverse momentum 
of the bachelor kaon, and their average, weighted by the number of events in each subsample, 
is computed. The difference from the result obtained in the default fit is 0.018 x 10~ 2 
(0.08 x 10~ 2 ), which is assigned as a systematic uncertainty. Finally, the CP measurements 
are repeated applying a geometrical requirement [261 t-o suppress possible detector edge 
effects. The associated systematic uncertainty is 0.13 x 10 -2 (0.05 x 10 -2 ). 

Using the known BC —)■ .J/i[) K ± CP asymmetry, A CP (B ± —> J/ipI^) = (+0.3 ± 
0.6) x 10~ 2 131, t-he asymmetries are measured to be A CP (B ± —s- rj'K ± ) = [—0.2 ± 
1.2 (stat) ± 0.1 (syst) ± 0.6 (norm)] x 10~ 2 and A CP (B ± ^ (jrK^) = [+1.7 ±1.1 (stat) ± 
0.2 (syst) ± 0.6 (norm)] x 10~ 2 , where the third uncertainty comes from the B±— > J/ijjK p 
CP asymmetry. These results are compatible with the hypothesis of CP symmetry and 
with the Standard Model predictions 12,27 


In conclusion, this Letter presents the first observation of the decay B ^ —$■ 7/7/, with a 
significance of 6.4 standard deviations, and the most precise measurements of CP-violating 
charge asymmetries in B ± —> ij'K ± and B ± —> (j)K ± decays. The latter result supersedes 
the previous LHCb measurement 1261. The measured B® —)■ 7/7/ branching fraction, 
B(B ° —>• rj'rj') = [3.31 ± 0.64 (stat) ± 0.28 (syst) ± 0.12 (norm)] x 10 -5 , agrees with the 
theoretical predictions. This newly observed B° s decay channel to a charmless CP eigenstate 
opens possibilities for further constraining the Standard Model with time-dependent CP 
asymmetry measurements. 
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